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We present experimental study of a topological excitation, interface, in a vertically vibrated layer of 
granular material. We show that these interfaces, separating regions of granular material oscillation 
with opposite phases, can be shifted and controlled by a very small amount of an additional sub- 
harmonic signal, mixed with the harmonic driving signal. The speed and the direction of interface 
motion depends sensitively on the phase and the amplitude of the subharmonic driving. 

PACS: 47.54.+r, 47.35.+i,46.10.+z,83.70.Fn 



Despite their ubiquity and many practical applications, 
an understanding of the fundamental dynamical behavior 
of granular materials remains a serious challenge |ij . One 
of the main obstacles for the development of a continuous 
description of granular flow is the difficulty in perform- 
ing quantitative experiments under controlled conditions. 
Testing of theoretical models of the basic excitations of 
the system is especially important. It has been shown 
recently (|-(^] that thin layers of granular materials sub- 
jected to vertical vibration exhibit a diversity of patterns 
which may play the role of such fundamental excitations. 
The particular pattern is determined by the interplay be- 
tween driving frequency / and the acceleration amplitude 
r = 4ir 2 Af 2 /g of the cell, where A is the amplitude of 
oscillation and g is the acceleration due to gravity. Pe- 
riodic patterns, such as squares and stripes, or localized 
oscillons vibrating with frequency f/2 appear at T w 2.4 
|^|-^|. At higher acceleration (r > 3.72), stripes and 
squares become unstable, and hexagons appear instead. 
Further increase of acceleration replaces hexagons with a 
non-periodic structure of interfaces separating large do- 
mains of flat layers oscillating with opposite phase with 
frequency f/2. These interfaces were called kinks in 
Ref. @||||]. These interfaces arc cither smooth or " deco- 
rated" by periodic undulations depending on parameters 
j|,U. For r > 5.7, various quarter-harmonic patterns 
emerge. Several theoretical approaches including molecu- 
lar dynamics simulations, order parameter equations and 
hydrodynamic-type models, have been proposed recently 
to describe this phenomenology, see e.g. [MSI. 

In this Letter we present an experimental study of the 
dynamics of interfaces in a vibrated thin layer of gran- 
ular material. We find that an additional subharmonic 
driving results in a controlled displacement of the inter- 
face. In the absence of subharmonic driving, the interface 
drifts toward the middle of the cell. When the subhar- 
monic frequency f/2 is slightly detuned, the interface 
moves periodically about the middle of the cell. The 
present experimental results are in agreement with theo- 



retical predictions M. 

Interfaces in a granular layer separate regions of gran- 
ular material oscillating with opposite phases with re- 
spect to the bottom plate of the vibrating cell. These two 
phases are related to the period-doubling character of the 
flat layer motion at large plate acceleration. Since an in- 
terface separates two stable symmetric dynamic phases, 
it can be interpreted as a topological defect, similar to 
a domain wall in ferromagnets separating regions of op- 
posite magnetization ]To| . Interfaces can only disappear 
at the walls of the cell or annihilate with other inter- 
faces. The existence of the interfaces can be understood 
from the following consideration. Since grains in a layer 
lose their kinetic energy in multiple inter-collisions dur- 
ing landing at the plate, the behavior of a granular layer 
can be compared with the dynamics of a fully inelastic 
ball bouncing on a plate vibrating with amplitude A and 
frequency / [p|JTl]]. In this case, for driving with accel- 
eration r less than To ss 3.72 the ball lifts to the same 
height at each cycle |3|,[l2| ■ Above Tq the motion exhibits 
period-doubling, i.e. the heights of elevation alternate at 
each cycle. As a result, for T > To, the two states of 
the bouncing ball, differing by the initial phase, would 
coexist. If the analogous states of the bouncing flat layer 
co-exist in different parts of the cell, they have to be 
separated by an interface. These interfaces, found ex- 
perimentally ||[||| and theoretically H , are flat for high 
frequency drives and show transverse instability leading 
to periodic decoration at lower frequencies (see Fig. 1). 
In an infinite system a straight interface is immobile due 
to the symmetry between alternating states: the motion 
of flat layers on both sides of the interface is identical with 
a phase shift tt. Additional driving at the subharmonic 
frequency f± = f/2 will break the symmetry between 
domains with opposite phases. Depending on the phase 
of the additional driving, with respect to the phase of 
the primary driving, the relative velocity of the layer and 
the plate at collision will differ on different sides of the 
interface, and the material on one side will be lifted to 
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a larger height than on the other side. As a result, the 
direction of the interface motion can be controlled by the 
phase <&. The speed at a given $ is determined by the 
amplitude of the subharmonic acceleration 7. 

In our previous work we have developed a phcnomeno- 
logical description of the pattern formation in thin layers 
of granular material [^j^] . On the basis of our order pa- 
rameter model we have predicted JoJ that for small values 
of 7, the interface velocity V = dX/dt, X being the po- 
sition of the interface, is a periodic function of $, 

^ =F sin($-$ ) (1) 

where Vo — cry, and susceptibility a and phase shift $0 
are functions of the driving amplitude A and frequency 
/ which can be estimated using the model Ref. || . 

We performed experiments with a thin layer of gran- 
ular material subjected to two-frequency driving. Our 
experimental setup is similar to that of We used 

0.15mm diameter bronze balls, and the layer thickness in 
most of our experiments was 10 particles. We performed 
experiments in a circular cell of 15.3 cm diameter and 
in a rectangular cell of 4 x 12 cm 2 , which were evacu- 
ated to 2 mTorr. We varied the acceleration T and the 
frequency / of the primary driving signal as well as accel- 
eration 7, frequency /1 and the phase $ of the secondary 
(additional) driving signal. The interface position and 
vertical accelerations were acquired using a CCD camera 
and accelerometers and further analyzed on a Pentium 
computer. The lock-in technique was used to measure 
the accelerations at / and //2 frequencies in real time. 

In the absence of additional driving the interface drifts 
toward the middle of the cell (see Fig. 1). We attribute 
this effect to the feedback between the oscillating gran- 
ular layer and the plate vibrations due to the finite ra- 
tio of the mass of the granular material to the mass of 
the vibrating plate. Even in the absence of subharmonic 
drive, the vibrating cell can acquire a subharmonic mo- 
tion from the periodic impacts of the granular layer on 
the bottom plate at half the driving frequency. If the 
interface is located in the middle of the cell, the masses 
of material on both sides of the interface are equal, and 
due to the anti-phase character of the layer motion on 
both sides, an additional subharmonic driving force is 
not produced. The displacement of the interface X from 
the center of cell leads to a mass difference Am on op- 
posite sides of the interface which in turn causes an ad- 
ditional subharmonic driving proportional to Am. In a 
rectangular cell Am oc X. Our experiments show that 
the interface moves in such a way to decrease the subhar- 
monic response, and the feedback provides an additional 
term —X/t in the r.h.s. of Eq. ([!]), yielding 

^ = -X/t + V sin($ - $„) (2) 

The relaxation time constant r depends on the mass ratio 
(this also holds for the circular cell if X is small compared 



to the cell radius). Thus, in the absence of an additional 
subharmonic drive (7 = 0), the interface will eventually 
divide the cell into two regions of equal area (see Fig. 1). 

In order to verify this model we performed the follow- 
ing experiment. We positioned the interface off center 
by applying an additional subharmonic drive. Then we 
turned off the subharmonic drive and immediately mea- 
sured the amplitude /1 of the plate acceleration at the 
subharmonic frequency p3| . The results are presented 
in Fig. 2. The subharmonic acceleration of the cell de- 
creases exponentially as the interface propagates to the 
center of the cell. The measured relaxation time r of the 
subharmonic acceleration decreases with the mass ratio 
of the granular layer and of the cell with all other param- 
eters fixed. The mass of the granular layer was varied by 
using two different cell sizes while keeping the thickness 
of the layer unchanged. For the cells shown in Fig. 1, we 
found that the relaxation time r in the rectangular cell is 
about 4 times greater than for the circular cell (see Fig. 
2). This is consistent with the ratio of the total masses 
of granular material (52 grams in rectangular cell and 
198 grams in circular one). In a separate experiment the 
mass of the cell was changed by attaching an additional 
weight of 250 grams to the moving shaft, which weights 
2300 grams. This led to an increase of the correspond- 
ing relaxation times of 15-25 %. The relaxation time r 
increases rapidly with T (see Fig. 2, inset). 

When an additional subharmonic driving is applied, 
the interface is displaced from the middle of the cell. For 
small amplitude of the subharmonic driving 7, the sta- 
tionary interface position is X = VoTsin(<i> — $0), since 
the restoring force balances the external driving force. 
Fig. 3 shows the positions of the interface for various <E>. 
From such images we determined the equilibrium position 
X as function of <!> (Fig. 4a). The solid line depicts the si- 
nusoidal fit predicted by the theory. Because of the finite 
mass ratio effect, the amplitude of the measured plate ac- 
celeration fj, at frequency /1 also shows periodic behavior 
with $, (see Fig. 4b), enabling us to infer the inter- 
face displacement from the acceleration measurements. 
For even larger amplitude of subharmonic driving (more 
than 4-5 % of the primary driving) extended patterns 
(hexagons) re-emerge throughout the cell. 

The velocity Vo which the interface would have in an 
infinite system, can be found from the measurements of 
the relaxation time r and maximum displacement X m at 
a given amplitude of subharmonic acceleration 7, Vo = 
X m /r, see Eq.(||). We verified that in the rectangular 
cell the displacement X depends linearly on 7 almost up 
to values at which the interface disappears at the short 
side wall of the cell (see Fig. 5, inset). Figure 5 shows the 
susceptibility a = V0/7 as a function of the amplitude of 
the primary acceleration T. The susceptibility decreases 
with r. The cusp-like features in the T-dependence of a 
(and t, see inset to Fig. 2), are presumably related to 
the commensurability between the lateral size of the cell 
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and the wavelength of the interface decorations. 

We developed an alternative experimental technique 
which allowed us to measure simultaneously the relax- 
ation time r and the "asymptotic" velocity Vq- This was 
achieved by a small detuning Af of the additional fre- 
quency /i from the exact subharmonic frequency f/2, 
i.e. Af = fx — f/2 <C /. It is equivalent to the lin- 
ear increase of phase shift 4> with the rate 2ttA/. This 
linear growth of the phase results in a periodic mo- 
tion of the interface with frequency Af and amplitude 
X m = Vq/^t- 2 + (2ttA/) 2 (see Eq (§). The measure- 
ments of the "response functions" X m (Af) are presented 
in Fig. 6. From the dependence of X m on Af we can 
extract parameters Vq, a and r by a fit to the theoretical 
function. The measurements are in very good agreement 
with previous independent measurements of relaxation 
time r and susceptibility a. For comparison with the 
previous results, we indicate the values for r and a, ob- 
tained from the response function measurements of Figs. 
2 and 5 (stars). The measurements agree within 5 %. 




FIG. 1. Stationary position of the interface in circular, di- 
ameter 15.3 cm (a) and rectangular 4 x 12 cm (b) cells for 
sinusoidal driving force with / = 40 Hz and T = 4.1. Layer 
thickness 10 particle diameters. 

In summary, the position of a vertically vibrated gran- 
ular layer can be controlled by a very small accelera- 
tion applied at the subharmonic frequency (of the order 
of 0.1% of the primary harmonic acceleration). The di- 
rection and magnitude of the interface displacement de- 
pend sensitively on the relative phase of the subharmonic 
drive. Our measurements confirm the theoretical predic- 
tions made on the basis of the order parameter model 



Q . We observed that period-doubling motion of the flat 
layers produces subharmonic driving because of the finite 
ratio of the mass of the granular layer and the cell. This 
in turn leads to the restoring force driving the interface 
towards the middle of the cell. 
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FIG. 2. Amplitude of the subharmonic acceleration /j, of the 
cell averaged over 4 measurements vs time for the interface 
propagating to the center of the cell for T — 3.97 and / = 40 
Hz. Circles/squares correspond to the circular/rectangular 
cells, open/closed symbols correspond to light/heavy cells, re- 
spectively. Heavy cells differ from light ones by an additional 
weight of 250 g attached to the moving shaft. Solid lines show 
exponential fit /j, ~ exp(t/r) + const. Inset: r vs V for light 
rectangular cell. 
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FIG. 3. Equilibrium position of interface for $ = 80 (a); 
$ = 170° (b); <E> = 260° (c) for V = 4.1, / = 40Hz, /1 = f/2, 
and 7 = 0.6% of P in a rectangular cell. 
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FIG. 4. (a) Equilibrium position X (b) amplitude of mea- 
sured subharmonic acceleration n as functions of phase $. 
Circular cell, F = 4.1, / = 40 Hz, 7 = 1.25% of P. 
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FIG. 5. Susceptibility a = V0/7 vs T at / = 40Hz, rectan- 
gular cell. Inset: Displacement X as function of 7 at $ = 260°. 
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FIG. 6. Maximum displacement X m from center of rect- 
angular cell as function of frequency difference Af = /1 — f/2 
for / = 40 Hz, and for T = 3.97 (circles) and V — 4.1 (dia- 
monds). Dashed lines are fit to X m = V /v/T- 2 + (27rA/) 2 . 
The values of a and r obtained from the fit are also indicated 
in Figs. 2 and 5 (stars). 
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